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The present knowledge about hadrons produced in relativistic heavy ion colhsions is compatible 
with chemical freeze-out happening when the energy density divided by the particle density reaches 
the value of 1 GeV. This observation is used to determine the energy dependence of the chemical 
freeze-out parameters Tch and /i^ for beam energies varying between 1 and 200 A-GeV. The con- 

■ sequences of this energy dependence are studied for various particle ratios. Predictions for particle 
' ratios at beam energy 40 A GeV are presented. The conditions for thermal freeze-out are also de- 
, termined. These correspond either to an energy density of 45 MeV/fm"^ or to a particle density of 

^ ■ 0.05/fm^ 

^3 : 
^ , 

o ■ 

(N , 

CsJ ■ I. INTRODUCTION 

> : 

' It is the purpose of the present paper to determine the energy dependence of the chemical freeze-out parameters, 
, namely, the temperature Tch and the baryon chemical potential fi'j^ for beam energies in the range between 1 and 200 A- 

■ GeV 0]. This covers the full range of energies available at the GSI/SIS, BNL/AGS and the CERN/SPS experiments. 
These parameters determine the particle composition of the hadronic final state. Below this temperature inelastic 

0^ ' collisions between hadrons are no longer important and the hadronic abundances remain unchanged. This information 
can then be used to make predictions for beam energies that are in between those of the above experiments and also, 
. eventually, to extrapolate results to higher beam energies. 

To determine in a relatively unambiguous way the chemical freeze-out values in relativistic heavy ion collisions it 
is best to use ratios of integrated particle yields. Such ratios are not very sensitive to the dynamics of the underlying 
processes, one of the reasons being that integrated particle yields are Lorentz invariant, i.e. a boost in the transverse 
^ , direction affects the momentum distribution of particles but not their number. The restriction to integrated particle 
!• yields therefore minimizes the model dependence of the freeze-out parameters. This of course not true for particle 
.'^ yields restricted to narrow kinematic regions since these will be more dependent on the dynamics of the process. 
^ During the past few years such analyses have been made for the CERN/SPS [|-||, the BNL/AGS and the 

5^ GSI/SIS data. A discussion and review of the various results can be found in Rcf. [|6j. 

The momentum distribution of particles produced in relativistic heavy ion collisions does not show the dependence 
predicted by a statistical distribution. The longitudinal direction bears no resemblance to an exponential, and in 
the transverse direction substantial deviations from an exponential fit with a universal slope corresponding to the 
temperature exist. 

In section 2 we discuss the general case of hydrodynamic flow starting from the Cooper- Frye formula As a 

particular case we consider boost invariant longitudinal flow accompanied by flow in the transverse direction. 
We next discuss the case where the longitudinal flow is reproduced by a superposition of fireballs. In all these cases 
particle ratios can be calculated as if they were given by a static Boltzmann distribution. 

In section 3 we discuss the exact conservation of strangeness in statistical models. This is necessary when the 
number of strange particles is very small as is the case in e.g. the GSI/SIS energy range. 

In section 4 we use the phenomenological observation that the freeze-out parameters all correspond to an average 
energy per hadron of 1 GeV, independent of the manner in which the system was created, to extract the energy 
dependence of Tch and /i^''. This dependence is then used to track the behavior of various hadronic ratios as a 
function of energy and to present predictions for a beam energy of 40 GeV. We also determine the conditions for 
thermal freeze-out in the energy range between 1 A-GeV and 200 A-GeV. The values obtained previously ]l9|-p7[ are 
compatible with an energy density e « 45 MeV/ini^ or a particle density n « 0.05 fm^"^. It is at present not possible 
to distinguish between these two possibilities. 
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In section 5 we summarize our results. 



II. RATIOS OF THERMAL PARTICLE ABUNDANCES 

In the thermal model particle densities are determined by a statistical distribution, denoted where the index i 
refers to the type of hadron, e.g i = tt"*", , ... etc. In the Boltzmann approximation is given by 

' V (2^)3" 

= 5m2Tif2(m,/r)e^'/^ (2.1) 

This expression applies for the density of particles inside a fireball at rest having a temperature T and particle i has 
a chemical potential fii . If the system is in chemical equilibrium then the possible values of fj,i are determined by the 
overall baryon chemical potentials for baryon number, /.ig, strangeness, ^^'^ charge fiQ. The value of fig is fixed 
by giving the overall baryon number density ng, and that of fis by fixing the overall strangeness to zero, the value of 
fiQ is fixed by giving the neutron surplus. If the number of initial state neutrons is equal to the number of protons as 
is the case in e.g. S — S then fig = 0, for Pb — Pb hq is small and negative. This means that for given values of T 
and pLB one tunes the remaining parameters, /.ts, and /ig in such a way as to ensure strangeness neutrality and the 
correct isospin composition of the system. 

At chemical freeze-out the particle density of hadrons of a given type is determined by summing over all hadronic 
densities multiplied by the appropriate branching ratio, e.g. for tt"*" one has 

^^niBr{i ^ Ti+). (2.2) 



We have included in Eq. 2.2 the contributions from all particles and resonances with masses up to 2 GeV [ p8[ . 
Thermal freeze-out happens when the elastic collisions between hadrons cease. Since it is obvious that the Boltzmann 
distribution is not a good description of the momentum distribution of particles in relativistic heavy ion collisions we 
will now discuss modifications which lead to a more realistic description of the particle spectra. 

A. Cooper - Prye Formula 

In a hydrodynamic description taking into account flow in the longitudinal and transverse directions, the final 
state particles will leave the hadronic gas at freeze-out time. The momentum distribution of particles is given by the 



Cooper-Frye formula |17 



^^-(2?)^//^^'^^^'''^- ^'-'^ 

where the integration has to be performed over the freeze-out surface described by cr^. Its direction is perpendicular 
to the surface and its magnitude is determined by the size of the freeze-out surface. For the temperatures under 
consideration it is safe to neglect quantum statistics and we will therefore work with the Boltzmann distribution from 
now on; the generalization to Fermi-Dirac or Bose-Einstein statistics is straightforward. We thus have 

f{x,p)=eMi.-P-^ + ^^)|T], (2.4) 

where T, /i and u*^ are the (space-time dependent) temperature, chemical potential and four-velocity, respectively 
pgj. As an example, for a static fireball, the freeze-out surface is given by 



da^ = (d3a;,0), (2.5) 



so that we obtain the expected result 



After integrating the differential distribution over all momenta one obtains the total number of particles ]l|,|30t 
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for a Boltzmann distribution this becomes 



(2.7) 



(27r)3 J, 
And as a consequence : 



A7:Tm'K,[^)] (2.8) 



(2.9) 



Thus effects of hydrodynamic flow cancel out in the hadronic ratio, provided there are unique freeze-out temperature 
and chemical potential. 

B. Rapidity plateau and transverse expansion 

For boost-invariant cylindrical expansion along the z-axis the correct variable to use is the proper time r (recall 
df^x = TdTdyrdrd(j)), we have 

da^ = (Tdyrdrd(j),f rdrdyrdcj))^ (2-10) 

where the second component is in the f direction, i.e. perpendicular to the surface of the cylinder. For the case where 
the flow is azimuthally symmetric, i.e. when an average is made over all events or when only head-on collisions are 
considered, one has therefore l3l],p2l 



/ dN^ \ ^ g_ 
\dymtdmt J tt 



r dr Tp[r) 



PfSinhj/A f mt cosh yt 



Tch ) V T^ch 

5z),,(«|?i»i),.„(«^) 

where Tp(r) refers to the freeze-out time which in general depends on r, so that the center of the cylinder freezes out 
before the surface. 

After integration over the transverse mass tot this leads to 

— — ^ - \ r dr TF{r) 

dy ) y=Q TT J a 

cosh(yT) - sinh(2/T)l mjT^nK^ e^^'^^K (2.12) 



If the freeze-out temperature is everywhere the same on the freeze-out surface this leads to the surprisingly simple 
result that the ratios are unaffected by the hydrodynamic flow 



{dNj/dy)y^Q m]T,hK2 (f^) e^-^/T^^ 



-| (2^3) 



This is a case where the momentum distribution in no way resembles a Boltzmann distribution : the rapidity distri- 
bution is flat and the transverse momentum distribution is affected by transverse flow, yet the integrated ratios are 
the same as those of a static Boltzmann distribution. 
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C. Superposition of Fireballs. 



In order to reproduce more closely the observed rapidity distribution it is useful to consider a superposition of 
fireballs along the rapidity axis |Q. The resulting particle density (integrated over the particle rapidity y) is given 
by 

n,= dy dYpB p{YFB)-7^{y-YFB) (2.14) 
J-oo J-Y dy 

where Yfb is the position of the fireball and p is the distribution of fireballs and Y is the largest value of the position 
in rapidity space of a fireball. In the above integral it is possible to interchange the integration limits so that one 
obtains 



rii = n. 



l-Y 

In a ratio the integral over the distribution of fireballs cancels out so that one is left with 



j ^dYpB p{Yfb) (2.15) 



^ = i (2.16) 



With the conclusion that the ratio is the same as if it were given by a Boltzmann distribution. Note that the 
cancellation is only possible if all fireballs are of a similar nature, i.e. they should all have the same temperature. 



III. EXACT STRANGENESS CONSERVATION 



In relativistic heavy ion collisions the energy and the number of hadrons in the final state is large enough to justify 
the use of the grand canonical ensemble. However if the temperature is very low, as is the case for the data from 
GSI/SIS, then the number of strange particles in the final state is very small and it becomes necessary to take this 



into account by using exact strangeness conservation 1 34 - 4l|] . The grand canonical ensemble only enforces strangeness 
conservation on average and allows for fluctuations around strangeness zero. This is not good enough if the number 
of strange particles is very small. The same holds for applications of statistical descriptions to a small system like the 
ones created in p — p or e"*" — e~ collisions as has been done recently in Ref . . 

We therefore treat strangeness exactly, especially in the GSI/SIS energy range but continue to use a grand canonical 
description for baryon number and for charge since the number of baryons and charged particles is always large. 

To restrict the ensemble summation to a fixed value of the strangeness 5, one performs the following projection 

Zs = ^ fj d(p e-*^^ Z(T, As, As, Aq) (3.1) 



where the fugacity factor Ag has been replaced by 



As = e^^. (3.2) 



and Z is the standard grand canonical partition function. Negecting for the moment the contribution of multi-strange 
particles like, S's and il's one can write the partition function as 

Zs=o - ^ rf0 exp {TVo + N^e't' + iV^ie"^^} (3.3) 

here A'^i stands for the sum of all single particle partition functions with strangeness plus one: 

7V_i EEiVA + iV^- + ... (3.4) 
while A^_i stands for the sum of all single particle partition functions with strangeness minus one: 

Ni=N-^ + Nk + .- (3.5) 
As an illustration we quote the explicit form of A'a assuming Boltzmann statistics 
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J d^pexp{{-EA + ^iB)/T), 



(3.6) 



It is clear from the above equation that baryon number is being treated grand canonically. 
To calculate the partition function more explicitly we expand each term in a power series: 



Zs=o = Zo— y^V] — -iV™JV!!i exp(im0) exp(- 

^TT /n ^ — ^ ^ — ^ ml Til 



(3.7) 



where Zq is the standard partition function for all particles having zero strangeness. Performing the integration over 
(h we are left with 



oo ^ 



(3.8) 



where one recognizes the series expansion of the modified Bessel function 

Zs=o = ZoIo{xi) 



(3.9) 



where Xi = 2^JN\NZi. The particle densities can be deduced in the standard way from the partition function. As 
an example we quote the number of kaons 



K ~ 



The total energy density can be calculated from 



similarly the particle number is given by 



{E) = T'—\nZs=o 



(N) —\nZs=o 



(3.10) 



(3.11) 



(3.12) 



Ai = l 



Of special interest is the small volume limit (or, more correctly, the small particle number limit). In this case the 
partition function can be expanded in a power series in which case the integration over (j) gives zero for all case where 
the oscillating factors do not match, i.e. 



Zs=o « 1 



1 

2^ 



(27r)3 
d^p 
(27r)3 



E/T+iij, 



-E/T-i4> 



(3.13) 



Thus it is clear that one needs at least two terms before one gets a non- vanishing result. For small particle numbers 
this leads to the well-known suppression of strangeness. 



IV. CHEMICAL AND THERMAL FREEZE-OUT PARAMETERS. 



A. Chemical Freeze-Out 



The considerations presented in the previous two sections make it possible to extract in a reliable way the values 
of the chemical freeze-out parameters T^h and /i^' from ratios of integrated particle yields. This procedure has been 
followed repeatedly over the past few years and results obtained from CERN/SPS, BNL/AGS and GSI/SIS data are 
shown in Fig. 1 and in Table 1, we also show the results obtained in Ref. from e"*" — e~ and p — p collisions. A low 
energy beam creates a hadronic gas with a corresponding low value of T^h but a high value of the baryon chemical 
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potential. This reflects the fact that such a system is predominantly made up of nucleons and contains relatively few 
mesons. At the other end, the CERN/SPS beam creates a hadronic system which has a high Tch but a very low value 
of /i^, this is because the hadronic system created is very rich in mesons and the proportion of baryons diminishes 
accordingly. It has been noted in Ref. |44| that, despite the wide variation in beam energies, all these points have in 
common that the energy density divided by the total particle density is 1 GeV, independently of whether the system 
was created in Ni — Ni at 1 A-GeV or in S* — S* at 200 A-GeV. The energy density has been calculated using Eq. 
(3.11) and the particle density was calculated using Eq. (3.12). This means that, independently of how the system 
was created from the moment the average energy per hadron becomes less than 1 GeV, inelastic collisions cease to 
be important and the chemical composition of the final hadronic state is fixed. We consider this to be extremely 
remarkable. 

Since the beam energy is increased smoothly as one follows the universal freeze-out curve, Fig. 1, it is possible to 
determine the energy dependence of the chemical freeze-out parameters, Tch and this will be done explicitly in 
the next section. We first discuss some of the implications of the freeze-out curve. 



B. Non-Relativistic Nature of Hadronic Gas 



The phenomenological observation that the chemical freeze-out always corresponds to an average energy per hadron 
of 1 GeV per hadron is easy to understand in the GSI/SIS energy range. Since the low temperature makes the hadronic 
gas a non-relativistic one, nucleons dominate and only very few other particles are present, one can therfore use the 
approximation 

|||«m^. + ^T«lGeV (4.1) 

which immediately reproduces the value of 1 GeV per hadron. 

At higher energies the hadronic gas becomes more and more a mesonic in nature and pions and p mesons dominate. 
It is of interest to compute the average mass of particles that make up the gas. This mass decreases very slowly, 
starting from the nucleon mass to approximately the p-meson mass. It is shown explicitly in Fig. 2. All hadrons 
in the gas are thus to a good approximation non-relativistic, except for pions which dominate for low values of the 
baryon chemical potential. The non-relativistic character is further enhanced by looking at the value of the average 
energy along the freeze-out curve, this is shown in Fig. 2. Thus the following relation holds 

|||«(M) + ^T«lGeV. (4.2) 

The particle composition of the hadronic gas changes smoothly as one follows the universal freeze-out curve from low 
to high temperatures. At low temperatures the hadronic gas is dominated by nucleons as can be seen explicitly in 
Fig. 3. For high temperatures the hadronic gas is meson dominated, pions and rho mesons making up the largest 
fraction. 



C. Thermal Freeze-Out Parameters 



After chemical freeze-out, the particle composition inside the hadronic gas is fixed but elastic collisions still keep 
the system together until the final, thermal freeze-out. At this stage the momentum distribution of particles no longer 
changes and is then final. The transverse momentum spectra therefore determine the thermal freeze-out parameters. 



and results are shown in Fig. 4. We have to 
DStantial, this introduces an additional uncertainty 



This has been done by several groups in the past few years |l£ 
remember that systematic errors on the particle spectra can be sut 
on the values of the parameters. The baryon chemical potential at this stage cannot be determined directly. We 
assume that the transition from chemical freeze-out to thermal freeze-out proceeds in such a way that the entropy to 
baryon number remains conserved, ie. S/B = constant. These are indicated by dashed-dotted lines in Fig. 4. We 
note that all results from 1 A-GeV to 200 A-GeV correspond to a fixed energy density of 45 MeV/fm'^, i.e. if the 
energy density inside the hadronic gas drops below this value, then the system ceases to exist independently of how 
it was formed. Wc note that it is not possible at present to distinguish between a fixed energy density and a fixed 
particle density. Both curves are shown in Fig. 4. 
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V. DISCUSSION 



A. Energy Dependence of Chemical Preeze-Out Parameters 

An inspection of Fig. 1 suggests several possibilities for extracting the dependence of the thermal parameters on 
the beam energy. It would e.g. be possible to plot the values of Tch as a function of the beam energy and simply 
interpolate between them. We follow a less direct method which relies on the use of results contained in Fig. 5. 
The values shown in this figure have been taken from the work of Gazdzicki |4^:^^ which uses a compilation of 
experimental results on the pion multiplicity divided by the number of spectator nucleons, (tt) /Apart- As can be seen 
from Fig. 5, in the relevant energy region this ratio increases approximately linearly with the beam energy or, more 
precisely with the variable 

\/s - ^Sthr (5.1) 

where Sthr is the threshold energy for pion production given by 

Sthr = '^'mN + m^. (5.2) 

This approximate linear dependence provides a convenient parameterization for extracting the energy dependence. 
A polynomial interpolation was also used in order to test the sensitivity of the parameterization. To extract the 
dependence on beam energy we combine Fig. 1 and Fig. 5. This is shown for one particular value in Fig. 6 where 
the crossing point between the two curves determines the value of the chemical freeze-out parameters. The resulting 
energy dependence of the freeze-out temperature T^h is shown in Fig. 7, as one can see, the value, at first increases 
rapidly with beam energy and then converges towards a maximum value of approximately 160 MeV. We note that 
this value is very close to the expected phase transition to a quark-giuon plasma as indicated by results from lattice 
gauge theory which give an upper limit of about 170 MeV for the critical temperature |^. The energy dependence 
of the baryon chemical potential is shown in Fig. 8. In this case the value at first decreases rapidly and then tends 
towards zero as the beam energy is increased. 

B. Hadronic Ratios. 

The energy dependence obtained for Tch and in the previous sub-section can be used to track various particle 
ratios as the beam energy increases. Of particular interest are the ratios of strange to non-strange particles since it 
is widely believed that these show the largest deviation from chemical equilibrium. Several analyses have indicated 
the need for an additional parameter, 75, which measures the deviation from chemical equilibrium of strange hadrons 
The j-K^ ratio is shown in Fig. 9, it increases smoothly from the low to the higher energies and reaches a 
maximum value of about 0.2 and stays approximately constant beyond this. We would like to refer to these values as 
the thermal model values for the I-k^ ratio . Comparing these results with the experimental ones we see that the 
thermal values nicely track the observed ratio in the GSI and BNL energy range. Additional values of the j-K^ 
ratio measured at mid-rapidity have become available recently from BNL/AGS in the range 2 - 10.7 A-GeV |5^. All 
reported values are in good agreement with the thermal model values. However the CERN/SPS value measured in 
S — S collisions is clearly below the predicted one. It has been argued recently that this deviation is due to 

the formation of a quark-gluon plasma in the initial state of the relativistic heavy ion collision. 

Next we consider the K'^ / K~ ratio. This ratio is extremely sensitive to the beam energy because it varies by two 
orders of magnitude over the energy range under consideration and is therefore a test of the thermal model. It is 
however not sensitive to the presence of 7s since it cancels out in the ratio. The comparison is shown in Fig. 10 and 
the agreement is good everywhere p5[ |. 

C. Predictions 

The results presented in the previous sections substantially increase the predictive power of thermal models. In 
the past it was not possible to use these models to calculate hadronic ratios in a new energy domain. In Table 2 we 
list predictions of the model for a beam energy of 40 A-GeV and compare them to experimental results obtained at 
CERN using Pb~Pb at 158 A-GeV. The main deviation can be found in ratios involving anti-baryons. This is because 
the baryon chemical potential is predicted to be substantially larger at 40 GeV than it 158 GeV, this suppresses all 
anti-baryons by an order of magnitude. 
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VI. SUMMARY 



We have shown that all results obtained in relativistic heavy ion collisions for beam energies between 1 A-GeV 
and 200 A-GeV can be summarized in a very simple way: when the average energy per hadron drops below 1 GcV, 
chemical freeze-out happens, the particle composition of the final state is fixed, and inelastic collisions cease to be 
important. This happens for Ni — Ni collisions at 1 A-GeV, for Pb — Pb collisions at 158 A-GeV and for S — S 
collisions at 200 A-GeV. It is thus independent of the beam energy, of the beam particle and by implication, of the 
size of the hadronic volume. Particle ratios for a beam energy of 40 A-GeV have been given in Table 2. We consider 
these predictions as highly reliable. 

The next stage of the evolution of the hadronic system is thermal freeze-out. Here the elastic collisions cease and 
the momenta of the final state particles are fixed. We have indicated that all available results are consistent with the 
fact that this happens when the energy density drops below 45 MeV/fm'^ or when the particle density is less then 
0.05/fm^. It is not possible yet to distinguish between these two possibilities with the presently available data. 

ACKNOWLEDGMENTS 

We acknowledge stimulating discussions with P. Braun-Munzinger, B. Priman, M. Gazdzicki, U. Heinz, W. 
Norenberg, H. Oeschler, H. Satz and J. Stachel. 



8 



iiucl-th/990301C 



on sabbatical leave at the Fakultat fur Physik, Universitat Bi elefeld, D-33615 B ielefeld, Germany. 
For a review see e.g. U. Heinz, J. Phys. G 25, 263 (1999) and |nucl-tli/9907060 . 
P. Braun-Munzinger, J. Stachel, J. P. Wessels and N. Xu, Phys. Lett. 365B, 1 (1995). 
P. Braun-Munzinger and J. Stachel, Nucl. Ph ys. A638, 3c (199 8). 
P. Braun-Munzinger, I. Heppe and J. Stachel, 
J. Letessier and J. Rafelski, J. Phys. G 25, (19"99)7 
C. Spieles, H. Stocker and C. Greiner, Z. Phys. G 74, 319 (1998). 

F. Becattini, M. Gazdzicki, and J. SoUfrank, Eur. Phys. Jour. G5, 143 (1998). 
I.G. Bearden et al., Nucl. P hys. A638, 419c (1998). 

G. D. Yen, M.I. Gorenstein, |nucl-th/9808012 . 

P. Braun-Munzinger, J. Stachel, J. P. Wessels and N. Xu, Phys. Lett. 344B, 43 (1995). 
J. Gleymans, D. Elliott, R. L. Thews, and H. Satz, Z. Phys. G74, 319 (1997). 
A.D. Panagiotou, G. Mavromanolakis and J. Tzoulis, Phys. Rev. G 53, 1353 (1996). 
J. Letessier, J. Rafelski and A. Tounsi, Phys. Lett. B328, 499 (1994). 

Gleymans, D. Elliott, A. Keranen and E. Suhonen, Phys . Rev. G 57 319 (1998). 
Gleymans, H. Oeschler and K. Redlich, ^lucl-th/980903l| , Phys. Rev. G 59 1663 (1999). 
SoUfrank, J. of Physics G 23, (1997) 1903 and references therein. 
Gooper and G. Frye, Phys. Rev. D 10, 186 (1974). 



J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 
H. Appelshauser et al.,NA49 GoUaboration, Nucl. Phys. 
S. Ghapman an d J.R. Nix, Phys . Rev. G 54, 866 (1996): 



J.R. Nix et al., tiucl-th/9801045[ 
M.M. Aggarwal et al., ^iucl-e:i^01009|. 209 (1998). 



A638, 91c (1998); Eur. Phys. J. G2, 661 
J.R. Nix, Phys. Rev. G 58, 2303 (1998). 



(1998). 



LG. Bearden et al., Phys. Rev. Lett. 78, 2080 (1998). 

E. Schnedermann, J. SoUfrank and U. Heinz, P hys. Rev. G 48, 2 462 (1993). 



J. Alam, J. Gleymans, K. Redlich and H. Satz, |nucl-th/9707042 . 

E.V. Shuryak, Nucl. Phys. A638, 207c (1998). 

B. Kampfer, J. Phys. G 23, 2001 (1997); hep-ph/961236. 

G. Gaso et al. (Particle Data Group), Eur. Phys. Jour. G3, 1 (1998). 

See e.g. D.H. Rischke, in "Hadrons in dense matter and Hadrosynthesis" , pp. 21-70, Ed. J. Gleymans, H.B. Geyer and 

E.G. Scholtz, Springer, (1999). 

D.H. Rischke, private communication. 

P.V. Ruuskanen, in "Quark-Gluon Plasma", pp. 519, Ed. R.G. Hwa, World Scientific, Singapore (1992). 
J. P. Blaizot and J.-Y. QlUtrault, in "Quark-Gluon Plasma", pp. 393, Ed. R.G. Hwa, World Scientific, Singapore (1992). 
J. SoUfrank and U. Heinz, in "Quark-Gluon Plasma", Vol. 2, Ed. R.G. Hwa, World Scientific, Singapore (1995). 
Hagedorn and K. Redlich, Z.Physik, G27, (1985) 541. 

Zalewski, Acta Phys. Pol. 28, (1965) 933; E. M. Ilgenfritz, and J. Kripfganz, Nucl. Phys. B62, 141 (1973). 
Rafelski and M. Danos, Phys. Lett. B97, 279 (1980). 

Redlich and L. Turko, Z. Phys. G5, 201 (1980); L. Turko, Phys. Lett. B104, 361 (1981). 

MuUer, " The Physics of the Quark-Gluon Plasma", pp. 91-104, Lecture Notes in Physics 225, Springer- Verlag (1985). 
Derreth, W. Greiner, H.-Th. Elze and J. Rafelski, Phys. Rev. G31, 360 (1985). 



R. 
K. 
J. 
K. 
B. 
G. 



M. Gorenstein, V.K. Petrov and G.M. Zinovjev, Phys. Lett. B106, 327 (1981). 
J. Gleymans, K. Redlich and E. Suhonen, Z. Phys. G 51, 137 (1991). 

F. Becattini, Z. Phys. G 69, 485 (1996); F. Becattini, and U. Heinz, Z. Phys. G 76, 269 (1997). 

E. Andersen et al., Phys. Letters, B433, 209 (1998). 

J. Gleymans and K. Redlich, Phys. Rev. Letters, 81, 5284 (1998). 

M. Gazdzicki, Z. Phys. G 66, 659 (1995). 

M. Gazdzicki and D. Rohrich, Z. Phys. G 65, 215 (1995). 

E. Laermann, Nucl. Phys. A610, Ic (1996). 

J. Rafelski, J. Letessier and A. Tounsi, Acta Phys. Pol. 27, 1037 (1996). 

F. Becattini, J. SoUfrank, K. Redl ich and H. Satz, Nucl. Phys. A638, 399c (1998). 
M. Gazdzicki and M. Gorenstein, |nucLth/9803462 . 

L. Ahle et al., E866 and E917 GoUaborations, Nucl. Phys. A6 38, 57c (1998). 
H. Dobler, J. SoUfrank and U. Heinz, Phys. Lett. B(in press), |nucLth/9904018 . 
B. Tomasik, PhD Thesis, Regensburg and to be published. 



9 



[54] L. Ahlo ct al., E866 and E917 Collaborations, Nucl. Phys. A638, 57c (1998). 
[55] F. Wang et al., NA49 Collab., liep-ph/982218. 
[56] F. Becattini, J. Phys. G 25, 287 (1999). 



10 



VII. FIGURE CAPTIONS 



Fig. 1 Chemical freeze-out parameters Tch and fi'g obtained at LEP, CERN/SPS, BNL/AGS and GSI/SIS. Refer- 
ences to the points are given in table 1. The full line corresponds to an energy density over total particle density of 1 
GeV. 

Fig. 2 Average mass along the freeze-out curve of Fig. 1 as a function of the chemical freeze-out temperature. 

Fig. 3 Change of the particle composition of the hadronic gas along the freeze-out curve of Fig. 1 as a function of the 
chemical freeze- out tem,perature. Shown are the fractions of nucleons and of pions. 

Fig. 4 Therm,al freeze-out (dashed and dashed-dotted lines) and chemical freeze- out (solid line) cunes. The dotted 
lines connecting the chemical and thermal freeze- out curves correspond to a fixed S/B (entropy /bary on number) ratio. 
Only references to the thermal freeze-out points are indicated explicitly. See Fig. 1 for references to the chemical 
freeze-out points. 

Fig. 5 Average number of pions divided by the number of participating nucleons as a function of the beam energy. The 
full line corresponds to a linear fit, the dashed line is a polynomial interpolation. 

Fig. 6 Intersection between the chemical freeze- out curve and the line corresponding to a given ratio of pions to Apart- 
The intersection determines the values of T^^ and /z^ for the beam energy corresponding to the indicated 'k/ Apart 
ratio. References to the chemical freeze- out points are given in Fig. 1. 

Fig. 7 Variation of the temperature at chemical freeze-out as a function of the beam energy. 

Fig. 8 Variation of the baryon chemical potential at chemical freeze-out as a function of the beam energy. 

Fig. ^Dependence of /ir^ ratio on beam energy along the chemical freeze- out curve of Fig. 1. Experimental results 

are indicated. See discussion in text. 

Fig. 10 Dependence of /K~ ratio on beam energy along the freeze-out curve of Fig. 1. Experimental results are 
indicated. 
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FIG. 1. Chemical freeze-out parameters T^h and ffg obtained at LEP, CERN/SPS, BNL/AGS and GSI/SIS. References to 
the points are given in table 1. The full line corresponds to an energy density over total particle density of 1 GeV. 
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FIG. 4. Thermal freeze- out (dashed and dashed-dotted lines) and chemical freeze- out (solid line) curves. The lines connecting 
the chemical and thermal freeze-out curves correspond to a fixed S/B (entropy /bary on number) ratio. Only references to the 
thermal freeze-out points are indicated explicitly. See Fig. 1 for references to the chemical freeze-out points. 
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FIG. 5. Average number of pions divided by the number of participating nucleons as a function of the beam energy. The full 
line corresponds to a linear fit, the dashed line is a polynomial interpolation. 
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FIG. 6. Intersection between the chemical freeze-out curve and the line corresponding to a given ratio of pions to Apart- The 
intersection determines the values ofTch and /ig for the beam energy corresponding to the indicated n/ Apart ratio. References 
to the chemical freeze- out points are given in Fig. 1. 
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Energy 


T^h [MeV] 


Ms [MeV] 


Reference 


SPS 








Pb + Pb 158 A-GeV 


170 ± 11 


270 ± 24 


i 


Pb + Pb 158 A-GeV 


165 


235 


i 


Pb + Pb 158 A-GeV 


145 ±5 


205 ±5 


i 


S + S 200 A-GeV 


180.5 ± 10.9 


220.2 ± 18.0 




S + Ag 200 A-GeV 


178.9 ±8.1 


241.5 ± 14.5 




s 9nn A .PpV 


1 71 n 


1 Rfl 04 


IHI 


o o zuu iA.-^^ev 


1 fin 9 -1- ^ 


lOo.U It ^.u 




a -V J ZUU 


ICC n -1- Pi n 


1 n -1- n 
1 / o.u =t o.u 




AGS 








Si + Au 14.6 A-GeV 


127 ±8 


485 ± 70 


P 


m 4- All 14 R A.PpV 


1 qri J_ 1 n 




Ilia 


m -1- All 14 fi A -CpV 


118 dz 12 


KOO _|_ CO 


IPOII 




IIU ztz 


c/in -1- on 


llril 


SIS 








Ni-Ni 1.9 A-GeV 


70± 10 


720 ± 50 




Ni-Ni 0.8 A-GeV 


48± 10 


820 ± 10 


PI 


Au-Au 1.0 A-GeV 


49 ±3 


825 ±8 


PI 


Ni-Ni 1.0 A-GeV 


51 ±5 


822 ± 10 




Ni-Ni 1.8 A-GeV 


54 ±3 


806 ±8 


PI 



TABLE I. Chemical freeze-out temperature Tch and baryon chemical potential in various collisions. 
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Energy 


Ty. [MeV] 


[MeV] 


Reference 


SPS 








Pb + Pb 158 K-GeV 


120 ± 12 


300 ± 35 




Pb + Pb 158 k-GeV 


95.8 ± 3.5 


400 ± 24 




Pb + Pb 158 k-GeV 


100 


400 




AGS 








Si + Au 14.6 A-GeV 


93.4 ±4.4 


554 ± 35 




Si + Au 14.6 A-GeV 


93 


540 


PI 







TABLE II. Thermal freeze-out temperature Tj and baryon chemical potential /^^ in various collisions. ^The value of pi?^ 
was estimated using isentropic expansion starting from the chemical freeze-out points. 
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Thermal 
Model 
40 GeV 



exp. data 
158 GeV 



Ref. 



(P - P)/h- 
P/P 
A/K« 
K+/K- 
A/A 
H^/A 



0.382 
0.006 
0.66 
2.34 
0.014 
0.12 
0.24 
0.038 
0.14 



0.228± 
0.055± 
0.65± 
1.85± 
0.128± 
0.127± 
0.180± 
0.227± 
0.46± 



0.029 

0.010 

0.11 

0.09 

0.012 

0.011 

0.039 

0.033 

0.15 



43 



TABLE III. Model Predictions for T = 
experimental results obtained at 158 GeV. 



140 MeV, fJ,B = 380 MeV corresponding to beam energy of 40 GeV compared to 
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